An environmental concern with hydraulic fracturing for shale gas is the risk of groundwater and surface water contamination. Assessing this risk partly involves the identification and understanding of groundwater–surface water interactions because potentially contaminating fluids could move from one water body to the other along hydraulic pathways. In this study, we use water quality data from a prospective shale gas basin to determine: if surface water sampling could identify groundwater compartmentalisation by low-permeability faults; and if surface waters interact with groundwater in underlying bedrock formations, thereby indicating hydraulic pathways. Variance analysis showed that bedrock geology was a significant factor influencing surface water quality, indicating regional-scale groundwater–surface water interactions despite the presence of an overlying region-wide layer of superficial deposits averaging 30–40 m thickness. We propose that surface waters interact with a weathered bedrock layer through the complex distribution of glaciofluvial sands and gravels. Principal component analysis showed that surface water compositions were constrained within groundwater end-member compositions. Surface water quality data showed no relationship with groundwater compartmentalisation known to be caused by a major basin fault. Therefore, there was no chemical evidence to suggest that deeper groundwater in this particular area of the prospective basin was reaching the surface in response to compartmentalisation. Consequently, in this case compartmentalisation does not appear to increase the risk of fracking-related contaminants reaching surface waters, although this may differ under different hydrogeological scenarios. 相似文献
The bonded discrete element model (DEM) is a numerical tool that is becoming widely used when studying fracturing, fragmentation, and failure of solids in various disciplines. However, its abilities to solve elastic problems are usually overlooked. In this work, the main features of the 2D bonded DEM which influence Poisson's ratio and Young's modulus, and accuracy when solving elastic boundary value problems, are investigated. Outputs of numerical simulations using the 2D bonded DEM, the finite element method, a hyper elasticity analysis, and the distinct lattice spring model (DLSM) are compared in the investigation. It is shown that a shear interaction (local) factor and a geometric (global) factor are two essential elements for the 2D bonded DEM to reproduce a full range of Poisson's ratios. It is also found that the 2D bonded DEM might be unable to reproduce the correct displacements for elastic boundary value problems when the represented Poisson's ratio is close to 0.5 or the long-range interaction is considered. In addition, an analytical relationship between the shear stiffness ratio and the Poisson's ratio, derived from a hyper elasticity analysis and applicable to discontinuum-based models, provides good agreement with outputs from the 2D bonded DEM and DLSM. Finally, it is shown that the selection of elastic parameters used the 2D bonded DEM has a significant effect on fracturing and fragment patterns of solids. 相似文献
This paper reports improvements to algorithms for the simulation of 3-D hydraulic fracturing with the Generalized Finite Element Method (GFEM). Three optimizations are presented and analyzed. First, an improved initial guess based on solving a 3-D elastic problem with the pressure from the previous step is shown to decrease the number of Newton iterations and increase robustness. Second, an improved methodology to find the time step that leads to fracture propagation is proposed and shown to decrease significantly the number of iterations. Third, reduced computational cost is observed by properly recycling the linear part of the coupled stiffness matrix. Two representative examples are used to analyze these improvements. Additionally, a methodology to include the leak-off term is presented and verified against asymptotic analytical solutions. Conservation of mass is shown to be well satisfied in all examples. 相似文献
Hydrogeochemical processes that would occur in polluted groundwater and aquifer system, may reduce the sensitivity of Sr isotope being the indicator of hydraulic fracturing flowback fluids(HFFF) in groundwater. In this paper, the Dameigou shale gas field in the northern Qaidam Basin was taken as the study area, where the hydrogeochemical processes affecting Sr isotope was analysed. Then, the model for Sr isotope in HFFF-polluted groundwater was constructed to assess the sensitivity of Sr isotope as HFFF indicator. The results show that the dissolution can release little Sr to polluted groundwater and cannot affect the εSr(the deviation of the 87 Sr/86 Sr ratio) of polluted groundwater. In the meantime, cation exchange can considerably affect Sr composition in the polluted groundwater. The Sr with low εSr is constantly released to groundwater from the solid phase of aquifer media by cation exchange with pollution of Quaternary groundwater by the HFFF and it accounts for 4.6% and 11.0% of Sr in polluted groundwater when the HFFF flux reaches 10% and 30% of the polluted groundwater, respectively. However, the Sr from cation exchange has limited impact on Sr isotope in polluted groundwater. Addition of Sr from cation exchange would only cause a 0.2% and 1.2% decrease in εSr of the polluted groundwater when the HFFF flux reaches 10% and 30% of the polluted groundwater, respectively. These results demonstrate that hydrogeochemical processes have little effect on the sensitivity of Sr isotope being the HFFF indicator in groundwater of the study area. For the scenario of groundwater pollution by HFFF, when the HFFF accounts for 5%(in volume percentage) of the polluted groundwater, the HFFF can result in detectable shifts of εSr(ΔεSr=0.86) in natural groundwater. Therefore, after consideration of hydrogeochemical processes occurred in aquifer with input of the HFFF, Sr isotope is still a sensitive indicator of the Quaternary groundwater pollution by the HFFF produced in the Dameigou shale of Qaidam Basin. 相似文献
In petroleum industry, the difference between pore pressure (Pp) and minimum horizontal stress Sh (termed the seal or retention capacity) is of major consideration because it is often assumed to represent how close a system is to hydraulic failure and thus the maximum hydrocarbon column height that can be maintained. While Sh and Pp are often considered to be independent parameters, several studies in the last decade have demonstrated that Sh and Pp are in fact coupled. However, the nature of this coupling relationship remains poorly understood. In this paper, we explore the influences of the spatial pore pressure distribution on Sh/Pp coupling and then on failure pressure predictions and trap integrity evaluation. With analytical models, we predict the fluid pressure sustainable within a reservoir before failure of its overpressured shale cover. We verify our analytical predictions with experiments involving analogue materials and fluids. We show that hydraulic fracturing and seal breach occur for fluid pressure greater than it would be expected from conventional retention capacity. This can be explained by the impact of the fluid overpressure field in the overburden and the pressure diffusion around the reservoir on the principal stresses. We calculate that supralithostatic pressure could locally be reached in overpressured covers. We also define the retention capacity of a cover (RC) surrounding a fluid source or reservoir as the difference between the failure pressure and the fluid overpressure prevailing in shale at the same depth. In response to a localized fluid pressure rise, we show that the retention capacity does not only depend on the pore fluid overpressure of the overburden but also on the tensile strength of the cover, its Poisson’s ratio, and the depth and width of the fluid source. 相似文献
For many decades most oil wells in Iran have produced using their natural flow potential and haven’t needed to be fractured. As time goes by, the reservoir pressure depletes and the need for hydraulic fracturing as a stimulation practice arises. Nonetheless there is no record of successful hydraulic fracturing in Iran.
The Bangestan reservoir with a suitable amount of oil in place and good rock reservoirs, has been selected for the present research work. In this work, the in situ stress profile was calculated by using the available petrophysical data. This is achieved by using poroelastic theory for the stresses, and the Mohr–Coulomb criterion to predict failure. The model leads to easily computed expressions for calculating the pressure required to maintain hydraulic fracturing. Then the appropriate depth for treatment was determined. The results indicate that Ilam and Sarvak formations could be good candidates for hydraulic fracturing. Then, for two layers, a hydraulic fracture was designed and the production was predicted and the Net Present Value (NPV) resulting from the fracture of both layers was investigated. 相似文献